Localized Random Lasing Modes and a New Path for Observing Localization 
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We demonstrate that a knowledge of the density-of-states and the eigenstates of a random system 
without gain, in conjunction with the frequency profile of the gain, can accurately predict the mode 
that will lase first. Its critical pumping rate can be also obtained. It is found that the shape of 
the wavefunction of the random system remains unchanged as gain is introduced. These results 
were obtained by the time-independent transfer matrix method and finite-difference-time-domain 
(FDTD) methods. They can be also analytically understood by generalizing the semi-classical 
Lamb theory of lasing in random systems. These findings provide a new path for observing the 
localization of light, such as looking for mobility edge and studying the localized states. 
PACS numbers: 42.25.Bs, 72.55.Jr, 72.15.Rn, 05.40.-a 
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Localization theory and laser theory were both devel- 
oped in the sixties. The propagation of quantum and 
classical waves in disordered media has been well under- 
stood , while the laser physics has been well estab- 
lished at the same time. It was always assumed 
that disorder was detrimental to lasing action. However, 
Letokhov [Q theoretically predicted the possibility of las- 
ing in a random system, called "random laser" . But only 
after the experimental observations by Lawandy et al. , 
the random laser systems have been studied intensively 
P-p^. Since then, many experiments were carried out 
that showed a drastic spectral narrowing pi and a nar- 
rowing of the coherent backscattering peak . Recently, 
new experiments [^,^ showed random laser action with 
sharp lasing peaks. To fully explain such an unusual be- 
havior of stimulated emission in random systems, many 
theoretical models were constructed. John and Pang ||lC| ] 
studied the random lasing system by combining the elec- 
tron number equations of energy levels with the diffusion 
equation. Berger et al. obtained the spectral and spa- 
tial evolution of emission from random lasers by using a 
Monte Carlo simulation. Very recently, Jiang and Souk- 
oulis |OLO|, by combining a FDTD method with laser 
theory ||2|,|3[, were able to study the interplay of localiza- 



tion and amplification. They obtained |T^,|T^ the field 
pattern and the spectral peaks of localized lasing modes 



inside the system. They were able to explain |12 13] the 



multi-peaks and the non-isotropic properties in the emis- 
sion spectra, seen experimentally ||,^ . Finally, the mode 
repulsion property which gives a saturation of the num- 
ber of lasing modes in a given random laser system was 
predicted. This prediction was checked experimentally 
by Cao et al. @ . 

One very interesting question that has not been an- 
swered by previous studies is what is the form of the 
wavefunction in a random laser system. How does the 
wavefunction in a random system change as one intro- 
duces gain? Does the wavefunction retain its shape in 
the presence of gain? Another very interesting point is if 
we can predict a priori, which mode will lase first? What 



will be its emission wavelength? If we understand these 
issues, we will be able to design random lasers with the 
desired emission wavelengths. In addition, we will be 
able to use the random laser as a new tool to study the 
localization properties of random systems. 

In this paper, we explore the evolution of the wave- 
function without and with gain, by the time-independent 
transfer matrix theory |15|-|l9||, as well as by time- 
dependent theory [^^. The emission spectra can be 
also obtained. In addition, we can also use the semi- 
classical theory of lasing [^|||] to obtain analytical results 
for the threshold of lasing, as well as which mode will 
lase first. This depends on the gain profile, as well as the 
quality factor Q of the modes before gain is introduced. 

Our system is essentially a one-dimensional simplifi- 
cation of the real experiments It consists of many 
dielectric layers of real dielectric constant £2 — 2.56eo (£0 
is the dielectric permeability of free space) of fixed thick- 
ness (&o = 100 nm), sandwiched between two surfaces, 
with the spacing between the dielectric layer assumed to 
be a random variable a„ = ao(l -I- W), where ag = 200 
nm and W has a random value in the range of [-0.7, 0.7]. 
We choose a 30-cell random system, as the first system 
of our numerical study. In Fig. la, we present the re- 
sults of the logarithm of the transmission coefficient as 
a function of frequency /. These results were obtained 
by using the transfer matrix techniques introduced in 
Ref. 15. Notice that we have three typical resonance 
peaks (denoted Pi, P2 and i-3) in the frequency range 
of 600 to 660 THz. As one can see from Fig. la, the 
linewidths of the three modes are different. P3 has the 
smallest linewidth and therefore the largest Q, while Pi 
has the largest linewidth and therefore the smallest Q. 
We have also numerically calculated the wavefunctions 
corresponding to these three peaks and indeed find out 
that the more localized wavefunction is the one with the 
larger Q. All the results above were obtained for the case 
without gain. 

According to the semiclassical theory of laser physics 
]^,^ , we generally use a polarization due to gain Pgam — 
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eox{^)E = eo{x'{i^) + ix"{^))E to introduce amplifying 
medium effects. Both x'(w) and x"i'-^) ^'^^ proportional 
to the outside pumping rate Pr and can be expressed by 
the parameters of the gain material |^ . 

To determine which peak will lase first, wc can again 
use the time-independent transfer matrix method (see 
Eq. (4) of Ref. 15) with a fvcqncncy-independent gain, 
which means the width of the gain profile is very large. It 
is well understood that time- independent theory |15-1^] 
for random lasers can be used to obtain the threshold for 
lasing. At threshold, the transmission coefficient T goes 
to infinity. In Fig. lb, we plot the logio(T) versus x" 
for the three peaks Pi, P2 and P3. In Fig. lb, the P3 
which has the largest Q has the smallest threshold for 
lasing. So, in the frequency-independent gain case, the 
transfer matrix method gives that the mode with largest 
Q will lase first. We have also used the transfer-matrix 
method with a frequency-dependent gain profile, given as 
a dotted line in Fig. la. We choose the central frequency 
as fa = oJaf^TT = 618.56 THz, which is exactly between 
the Pi and P2, and its width to be Afa = Aw(j/27r = 15 
THz. In this case, when we increase the pumping rate Pr, 
we find the P2 will lase first, then Pi. So two important 
conditions determine which mode will lase first, the first 
is the quality factor of the mode and the second is the 
gain profile. Experimentally, it is quite often that only 
part of the random medium has been pumped. Then a 
third factor comes in, i.e. the spatial overlap of the mode 
function and the gain region. 
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FIG. 1. (a) Logarithm of transmission coefficient vs fre- 
quency of a 30-cell random system with three typical peaks. 
The dotted line shows the frequency dependent of the gain 
profile. (b) Logarithm of the transmission coefficient vs 
x" for the three peaks shown in Fig. la with a fre- 
quency-independent gain, x" is proportional to the pumping 
rate Pr- 



The next issue we address is the shape of the wave 
function, as one introduces gain. In Fig. 2a, we present 
the amplitude of the electric field versus distance for 



/ = 626.6 THz, which corresponds to the peak P2 of 
Fig. la. In Fig. 2b, we show the wavefunction with 
near-threshold gain at the exact high-peak frequency of 
P2- Both Fig. 2a and 2b are obtained by the time- 
independent transfer matrix method. Notice that in the 
presence of gain, the shape of the wavefunction in Fig. 2b 
is almost same as that without gain in Fig. 2a. The only 
change is its amplitude which increases uniformly. Actu- 
ally, by keeping the incident amplitude the same, when 
we increase the gain from zero to the threshold value, 
we find that the amplitude of the wavefunction increases 
from a small to a very large value, but its shape remains 
almost same. 
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FIG. 2. Amplitude of the electric field of the 30-cell system 
vs x/xq, with So = 5 nm, for the peak P2 in Fig. la. (a) With- 
out gain, incident field Einc = IV/m; (b) With near-threshold 
gain ( Pr = 7.2 x lO^'/s ), incident field Ei^^ = IV/m; (c) The 
lasing field with over-threshold pumping rate (P,. = 2xl0'^/s . 
Both (a) and (b) are obtained by the transfer matrix method 
while (c) is obtained by the FDTD method and laser theory. 

To check if indeed this amazing property is also obeyed 
by the full time-dependent theory 



13 1 of the semiclas- 



sical laser theories with Maxwell-equations, we repeat our 
theory for this case too. While in the time-independent 
theory, every mode is independent and amplification is 
not saturated, this is not true for a real lasing system. In 
a real lasing system, modes will compete with each other. 
As was discussed in Ref. 12, in a short system, the first 
lasing mode will suppress all other modes, so we observe 
only one lasing mode even for a large pumping rate. This 
is exactly the case when we use the FDTD method to 
simulate our 30-cell system with an over-threshold gain 
whose gain profile is the same as given in Fig. la. At 
first, the electric field is a random one due to sponta- 
neous emission, then a strong lasing mode evolves from 
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the noisy background and a sharp peak appears in the 
emission spectrum after Fourier transform. The stable 
field profile of the FDTD calculation is given in Fig. 2c, 
where the wavefunction is quite the same as in Fig. 2a 
and 2b, but the amplitude is very large. The emission 
spectrum of this case is obtained, and gives a sharp peak 
very close to the resonant peak P2 of Fig. la. We also 
checked the wavefunctions as well as the lasing threshold 
as we shift the gain profile. If the central frequency of the 
gain profile is near P3, the P3 mode will lase first, and 
the shape of wavefunction of P3 remains unchanged. We 
have also checked the above ideas for a larger (L ^ ^, 
where ^ is the localization length) system. In such strong 
localized case, the form of the wavefunction remains un- 
changed as one introduces gain. 

The numerical results which are presented above can 
be explained by the semi-classical Lamb theory Q of laser 
physics. According to Lamb theory, the Maxwell equa- 
tion of the random laser system can be written as 



the mode, called pulling effect. For a well-defined mode, 
generally Qm ^ 1, we need a very small gain to lase. 
Then x'(a;,a;) ^ 1, the pulling effect is very weak, and 
Lu ~ r2,„, where f2,„ is the eigenfrequency of the mode. 
So the wavefunction of the lasing mode should be simi- 
lar to the eigenfunction of the mode. Theoretically, we 
can use the perturbation method to obtain lu and the 
wavefunction. 

Equation (3) is the time-dependent amplitude equa- 
tion. First we can use it to determines the threshold 
condition when —x"{uj) = 1/Qm- For our system, with 
homogeneous pumping, we have 
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4{UJ^-Qm.) 



(4) 



where Cq 



is a constant. Eq. (4) 



dE{x,t) 
dt 



, ,d^Eix,t) 92 
- l^oe(x) ^ — = -^0 



Where /io is the permeability of free space, E{x, t) is the 
electric field , and the dielectric constant e{x) is deter- 
mined by the random configuration of the system, a is 
not just the common conductivity loss, but can be inter- 
preted by the total mode loss including the surface loss 
of the system by radiation. The polarization due to gain 
Pgain{x^t) is the Same as the one defined above. 

We assume that the system satisfies the slow varying 
approximation (it is always satisfied if we care about the 
stable lasing state). So E{x,t) = Em{t)(j)m{x)exp{—iujt) 
(our later discussion shows that the separation of spa- 
tial and time parts of wavefunction is reasonable), where 
Emit) is the field amplitude, (pmix) is the normalized 
wavefunction, and uj is the frequency of the field. The 
surface loss of the mode is Gm — Sa^/Qm, where the 
Qm is the quality factor of the mode. Then, we get two 
equations H for the real and the imaginary terms of Eq. 
(1) 



indeed gives that the threshold value of lasing is inversely 
proportional to the quality factor Q. 
p . j-^e^ynd, Eq. (3) give us the stable amplitude of field 
ain is over threshold. Actually, the gain is 
saturable, x"(w) k AA^ oc 1/(1 -t- CalE^P) [|o), in real 
systems and in our FDTD model ||l^,|l3| . With an over- 
threshold gain, Em{t) will increase and the gain param- 
eter x"(w,£^m) will decrease until —x"iuJ,Em) = l/Qm, 
then the field is stable. So Eq. (3) also determines the 
amplitude of the stable field for over-threshold pumping 
cases. 



\/^4>mix) + Mo(e(a;) + eox'ix,uj))uj^4>7nix) = 



dEmjt) 
dt 



= -X"{u^) - 



eQUjE^jt) 
2e 



(2) 



(3) 



where e ~ J^^ e{x)dx/L is the spatial averaged di- 
electric constant inside the system, and x"(ti^) — 
Jq 4'mix)*x"ix,uj)(j)mix)dx/L IS the spatial averaged 
gain. The last integral is done to take into account of 
the overlap between the wavefunction and the spatial re- 
gion of gain. 

Equation (2) determines the field distribution, quality 
factor Qm, and the vibration frequency uj of the lasing 
mode. The term Sox'{x,U!) will cause the vibration fre- 
quency to shift away from the original eigenfrequency of 
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FIG. 3. Logarithm of the transmission coefficient T, the 
emission intensity 7e and the gain profile vs frequency / for a 
80-cell random systems. In the insert, wavefunctions of two 
modes are shown. 

Our numerical and analytical results clearly suggest 
that states of the random system with gain can easily 
lase, provided their Q-factor is large. These findings pro- 
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vide a new path for observing localization of light. Since 
localized states have large Q values, will lase with a small 
pumping rate. On the other hand, strongly fluctuating 
extended states have smaller Q values because of the ra- 
diation loss on the surface of the system and can lase only 
after a stronger pumping. In a real experiment, even in 
the presence of absorption, if the gain profile is close to 
the mobility edge, there is going to be a discontinuity 
in the critical pumping rate needed for lasing. Localized 
states will lase first at a low pumping rate, while extended 
states need a high pumping rate. In Fig. 3, we present 
the results of the DOS vs frequency / for a Id quasi- 
disordered system. On the insert two eigenfunctions are 
given, one localized at / = 468 THz and the one "ex- 
tended" at / = 485 THz. In addition the emission inten- 
sity If. vs / is also shown with gain, Pr — 10^ 1/s, which 
demonstrate that the localized modes will lase first. In a 
more realistic 3D case, we expect le vs frequency to have 
many peaks for the localized region and no peaks in the 
extended region for a give pumping rate It would 

be very interesting if this discontinuity can be observed 
experimentally. 

In summary, we have used the time-independent trans- 
fer matrix method and the FDTD method to show that 
all lasing modes are coming from the eigenstates of the 
random system. A knowledge of the eigenstates and the 
density-of-states of the random system, in conjunction 
with the frequency profile of the gain, can accurately 
predict the mode that will lase first, as well as its criti- 
cal pumping rate. Our detailed numerical results clearly 
demonstrate that the shape of the wavefunction remains 
unchanged as gain is introduced into the system. The role 
of the gain is to just increase uniformly the amplitude of 
the wavefunction without changing its shape. These re- 
sults can be understood by generalizing the semi-classical 
Lamb theory of lasing in random systems. These findings 
can help us unravel the conditions of observing the local- 
ization of light, as well as manufacturing random laser 
with specific emission wavelengths. 
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